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Abstract. For the two dimensional Ising model, we construct the adequate surface 
tension near criticality. The latter quantity has been shown [9] to play a central role 
in the study of phase coexistence in a joint limit where the temperature approaches 
the critical point from below and simultaneously the size of the system increases fast 
enough. 



1. Introduction 

The present paper is devoted to a fundamental quantity related to the phe- 
nomenon of phase coexistence, namely the surface tension. In particular, we are 
interested in its behavior when we approach the critical point. 
There is a general belief that the surface tension and other quantities become 
isotropic near the critical point. Indeed, a wide family of statistical physics models 
defined on discrete lattices and taken at the critical point are conjectured to con- 
verge to a conformally invariant measure when the lattice mesh size goes to zero. 
In the particular case of independent site percolation on the triangular lattice, this 
conjecture has been made rigorous in [25] and [4]. There, the continuous confor- 
mally invariant measure is described by the Schramm-Loewner Evolution [22]. 
Critical phenomena also appear when one is not strictly at the critical point. When 
considering the model " near criticality" , that is when one takes the thermodynam- 
ical limit and simultaneously approaches the critical point, then the asymptotic 
behavior of relevant quantities is also influenced by critical phenomena. Among 
these joint limits, there is a special one, sometimes called the "scaling limit", at the 
threshold between two different behaviors. 

On one side, if the temperature approaches the critical point fast enough, the 
model behaves as if the temperature is exactly at the critical point. In the case of 
independent percolation, such a behavior has been proved by Kesten in [17]. 
On the other side, when the temperature approaches the critical point slowly enough 
then the influence of criticality is different. In this case, the behavior of the model 
can be described by non-critical phenomena that are altered by criticality. Indeed 
in [9] , it has been proved that the 2d-Ising model taken at sub-critical temperatures 
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at the vicinity of the critical point and in boxes that are large enough, still exhibits 
phase coexistence. In such regimes, the Wulff crystal persists but criticality washes- 
out the anisotropy inherited from the geometry of the lattice and reduces the Wulff 
shape to an ordinary circle. 

The paper [9] required the construction of an adequate joint-limit surface tension. 
The proof of this construction has been sketched in [19], using heavily the explicit 
computations techniques of [11]. In the present paper, we provide a construction 
of the surface tension near criticality without adding an extra layer of explicit 
computations. Our results are, for the time being, restricted to the two dimensional 
Ising model on the square lattice. Indeed, the starting point of our analysis is the 
beautiful exact formula for the fixed temperature surface tension. 

Using the duality property of Kramers- Wannier [18], the surface tension of the 
2d-Ising model at an inverse temperature (3 > (3 C can be defined from the asymptotic 
behavior of the two point function at the dual inverse temperature (3 < (3 C [21]. To 
be more precise, for each (3 < f3 c , let us denote by nr the unique infinite volume 
Ising measure on the spin configurations a e { — 1, 1} Z . Then for each x € Z 2 and 
at each (3 > (3 C , the surface tension is the function defined by the following limit 

(1) Tf s{x) = - lim - log t/[cr(0)cT(nx)}, 

n — >oo fi 

where (3 < (3 C is related to (3 > [3 C by the duality relation 

sinh(2/3) sinh(2/3) = 1. 

It is well known [21] that the function Tg can be extended continuously into a norm 
onM 2 . 

An important particularity of the 2d-Ising model is its relation to the dimcr model 
[16]. Indeed, Kastelcyn discovered that the partition function of the Ising model 
can be represented as the generating function of a dimer model. This permitted to 
give an explicit formula for the partition function. 

Later, a judicious application of Kasteleyn's representation and a tricky asymptotic 
analysis of Tocplitz matrices enabled Mc Coy and Wu [11] to derive the precise 
asymptotic behavior of the two point function between the spins of two distant 
sites. From these computations, it is possible to give via (1) a beautiful formula for 
the surface tension. Their result states that for all (3 > f3 c , and x — (x\,X2) € Z 2 : 

(2) T p( x ) — xiarcsinh(y / 1 + s 2 x\) + a^arcsinh^ 1 + s 2 ^), 
where s solves 

yji + s 2 *! + y/l + * 2 4 = smh(2/?) + — L^. 

This formula is the only input from explicit computations that we use in this paper. 
In [11], the formulas describing the asymptotic behavior of the correlation function 
and the surface tension are not given in the form above. In the last section of 
this paper, we will show how the results of Mc Coy and Wu can be written in the 
form (2). As we will see, this rewriting reveals a simple connection between the 
geometry of the 2d-Ising Wulff shape and the large deviation rate function of the 
simple random walk on Z 2 . 

The formula (2) is the result of the limit (1), i.e., n — > oo but at a fixed (3. The 
main purpose of our work is to derive the asymptotic behavior of the quantity in 
(1) in a joint limit (3 1 (3 C and n — > oo. Let us state the main result of this paper. 
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1.1. The main result. 

Theorem 1. For all x G 1? and for any double sequence n | co and [3 I [3 C 
satisfying 

n 

(3) : (p-(3 c )^oo, 



logn 



we have that 



(4) -lim 1 log^[a(0)a(nx)]=4\x\, 

n,0 (fj - (3 c )n 

where \x\ is the ordinary Euclidean norm of x. 

In words, the above result states that in a joint limit satisfying (3), it is still possible 
to define a rescaled surface tension and that the norm associated to this surface 
tension is the ordinary Euclidean norm. The latter fact is an indication that the 
model studied in the regime (3) is asymptotically rotation-invariant. The regime 
of Theorem 1 is, up to a logarithmic correction, sharp. Indeed, in [2, 3, 26] it has 
been shown that in a regime where n((3 — f3 c ) = t stays constant, one gets that 

. i ,. u^\a(0)a(nei)] , . 

(5) lim 1/4 =q* ' 



where ei is the unit vector in the horizontal direction and rj(t) satisfies a Painleve 
equation. Thus, Theorem 1 delimits the threshold between the constant temper- 
ature case and the regime considered in (5). The existence of the joint limit (4) 
is proved in two steps. First, we use (2) to compute the limit (4) when we first 
take the limit n | oo and then the limit (3 [ (3 C . In the second step, which is 
our main contribution, we provide a probabilistic argument that shows that in any 
joint regime satisfying (3) the joint limit of (4) is well defined and equals the limit 
obtained in the first step. First, we give the easy step of the proof of Theorem 1. 

Proposition 2. Uniformly over x <G K 2 , we have that 

lim lim J° g ^(0ML^])] = Um Mx) = 

m^iToo (f3 - [3 c )n\x\ 0lf3c {[3-[3 c )\x\ 

where for x — (xi,:^), \nx\ — ([nxij, [^2]) and for i € {1,2} [nxi\ is the largest 
integer which is smaller than or equal to nxi. 

The proof of Proposition 2 is obtained directly from the formula (2) by taking the 
limit (3 I [3 C . 

1.2. Comparison with independent site percolation. Recently, a lot of 
progress has been made in the study of criticality in the context of planar in- 
dependent site percolation on the triangular lattice T. In this section, we compare 
the known results on the Ising model with their analogue in independent site per- 
colation. The percolation analogue of rp is given by 

Vp > 1/2 Vx e T tJx) = - lim - log^O <-> nx), 
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where p = 1 — p < 1/2 and where P p is the probability measure corresponding to 
the site percolation of parameter p on T. 

Proposition 2 is actually a very strong statement that implies the existence of the 
correlation length exponent v = 1 in a rather strong from. It shows that the 
correlation length, if correctly rescaled becomes isotropic near the critical point. 
This result is at present time impossible to obtain without explicit computations. 
In percolation, Smirnov and Werner [25] used the convergence of critical percolation 
to the Schramm Loewner Evolution [24, 4] and the scaling relations of Kesten [17] 
to prove that 

VxGT lim , l0 r p( fL =4/3- 
pll/2 Iog(p - 1/2) 

And the question whether the limit lim p ^i/ 2 (p — 1/2)~ 4//3 t p (x) exists is still open. 
The appearance in (5) of a function r](t) which is related to a Painleve equation is 
very striking. Such asymptotics are unavailable for percolation. However, an ade- 
quate modification of the critical scaling limit has been proposed in [5] to investigate 
a regime analogous to the one considered in (5). 

1.3. Organization of the paper. In the next section we introduce the necessary 
notations and prove Theorem 1. The third section contains a rewriting of the results 
(2) of Mc Coy and Wu that permit to make the link with random walk and to draw 
some heuristics that explain the appearance of isotropy near the critical point. 

2. Proof of the Theorem 

The proof relies basically on sub-additivity. The easier part is the upper bound: 

Lemma 3. For all i£Z 2 and for all double sequences n f oo and (3 J. (3 C , we have 
that 

limsup — — — log^ /3 [cr(0)(T(nx)] < — 4|a;|. 

Proof. Let us fix (3 < (3 C and i€Z 2 . By the FKG- inequality, the sequence 

(log/[a(0)a(nx)],n>l), 

is super-additive and thus, for every fixed n > and (3 < (3 C , 

(6) pP[a{Q)a{nx)\ < exp ( lim — log /J* [a (0) a (mnx)] ] = cxp (— utr{x)) . 

\m]oo m J 

The result follows from Proposition 2 that guarantees that 

lim 1 tr(x) = 4|d. 
WA 13- P c 

n 

To prove the lower bound, we consider a sub-additive quantity that approximates 
adequately the two point function. To construct this quantity, we introduce for 
each (3 > (3 C , the unit ball of the norm defined by Tp: 

U? = {yeR 2 : T P {y)<l}. 
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U 13 for increasing /?, the number represents the radius of the dashed circle. 



Next, for fixed [3 > f3 c , n > and i e Z 2 , we define the tangent line w^(nx) of 
Tfj{nx)U^ at the point nx. We will also denote by w^' + (nx) the half plane delimited 
by w^(nx) and which does not contain the origin. We will also need the following 
discrete approximation of w@ (nx) : 

(7) W (nx) ={ye w^' + (nx) C\I? : < (y - nx) ■ v p (x) < 1}, 

where v^(x) is the unit length vector which is perpendicular to w^(nx) and which 
points towards the interior of w^' + (nx). 
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The quantity that permits us to establish the lower bound is 
(8) g(nx)= J2 A^My)]. 

The study of g(nx) started in [10] and in [1]. The results of the last two cited papers 
where then refined [15, 6, 7, 8] to obtain Ornstein-Zernike asymptotics for the two 
point function. Even though we establish rougher estimates than Ornstein-Zernike 
asymptotics, the quantity g(nx) turns out to be important when analyzing the two 
point function near criticality. 

Lemma 4. For each fixed (5 > [3 C and ieZ 2 , we have that: 
i) The following limit exists 

lim--log Yl A^)a(y)}, 
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and is equal to the surface tension Tp{x). 
ii) There exists a positive constant K which does not depend on x,n nor in [3 such 
that 

Vn > 1 E /^k(OMy)] > exp(-nr /3 (x) - K). 

yeWP{nx) 

Proof. The proof is an application of Simon's correlation inequality [23] which states 
that for any two sites x,y £ Z 2 and for any set W separating the two sites x, y, the 
following holds 

(9) A°{x)<v)] < E A^M^Hz)a(y)]. 

z£W 

Note that for independent percolation a similar inequality holds thanks to the van 
den Berg-Kesten inequality. 
Fix n, m > and consider 

g((n + m)x)= ]T f/[a(0)a(y)}. 

y£WP({n+m)x) 

For every y e W^^n + m)x), the set W^fax) separates the site from y, so that 
by (9) we obtain 

g((n + m )x)< E E ^HO)<r(z)]Aa(z)a(y)}. 

yeW {(n+m)x) zeW' 3 (nx) 

By translation invariance of the measure we get 

g((n + m)x)< E E A^0Mz)}^[a(0)a(y - z)]. 

zeWf 3 (nx) yeWI 3 ((n+m)x) 

By (7), for every z <G WP(nx) and y e W^((n + m)x), we have that 

< (y — z — mi) • v l3 (x) < 2. 

Therefore, there exists a site s £ W"(mx) such that \s — (y — z)\ < 2. Using the the 
FK-rcprcsentation Theorem 1.16 of [13] and the finite energy property Theorem 
4.17 of [13], we get that 

/K0)a(y - z)] < eVMOM*)], 
where if is a positive constant independent of [3. Hence 

g{{n + m)x)<e K E E A^0)a(z)]^[a(0)a(s)} 

zeWf(nx) s£WP(mx) 

<e K g(nx)g{mx). 
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By standard sub-additivity arguments, we get that the limit 

lim - log g(nx) = -gp(x) 

njoo n 

exists and moreover, for every n > 1 

g(nx) > exp(-ng (x) - K). 
So it remains to prove that gp (x) — Tp (x) . We write 

J2 ^[*(OMy)]=A<r(PMnx)]+ J2 /WOMy)] 

yeWP(nx) yeWfi(nx)\{nx} 

Using (6) and the definition of ^(ra), we have 

J2 f/[a(0)a(y)}< ]T exp^r^y)) 

y€WP(nx)\{nx} y£W& (nx)\{nx} 



oo 



<J2 J2 exp(-r /3 (2/)). 



k — y£W0(7ix) 

k<T l3 (y)-T l3 (nx)<k + l 

To evaluate the last sum, we have to establish an upper bound for the cardinality 
of the set 

{y <E W (nx) : k < T fj (y) - Tp(nx) < k + l} , 

where k > 0. For this, we note that there exists two positive constants ci and Ci 
such that 

Vj/GZ 2 Cl \y\ < r f3 (y) < c 2 \y\. 

From there we conclude that 

{y e W^inx) : k < Tp(y) - T f3 (nx) < k + l} C 

c{y e W l \nx) : \y\ e ( % \k + Tp(nx)), c^(k + 1 + Tp{nx)) } } . 

Using the fact that Tp{nx) < cin\x\ we can bound the width of the interval con- 
taining \y\: 

a 1 1 (k + 1 + T{j(nx)) - c^ 1 (k + rp(nx)) < — — — (fc + c 2 n\x\) + — . 

Therefore, for n large enough we can find a positive constant C such that 
\{y e W f \nx) : k < Tp{y) - Tp(nx) < k + l}| < C(k + 1 + n). 

Hence 

oo 

V?[<r(0Mv)] <Cexp(-nT /3 (x)) ^(fc + 1 + n) e - fe 

yeW 1 (nx)\{nx} fc=0 

- c \— + T^W) cM - nTp{x)) - 

Thus 

lim - log V fiP[a(0)a(y)] = lim - \ogfjP[a{Q)a{nx)] = -t {x). 

n|oo n n|oo n 

□ 

It is crucial to note that in the case x = ae i: i e {1,2} where a e N\{0} and (ei, e 2 ) 
is a orthonormal basis of R 2 , then one can derive a slightly stronger result. 
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Lemma 5. For each fixed a E N \ {0}, (3 > [3 C and x = aei,i = 1,2, we have that: 
i) The following limit exists 



iim--io g J2 /KOMy)], 



and is equal to the surface tension Tfj{x). 
ii) For all n > 

(J?[o-(0)o-(y)} > expi-nrpix)). 

y£w@ (nx) 

Proof. The proof is almost the same than the proof of Lemma 4, except that 
W^(nx) can be replaced by the line w^(nx) then by translation invariance we have 
that for every integers n and m and for each z € w^{nx) and y £ w^((m + n)x), 
we have that i/-ze w^{mx). This turns (g(nx), n > 1) into an exact sub- additive 
sequence. This property permits to prove the enhancement ii) for such a particular 
x. □ 

The last two lemma indicate that (8) is a nice quantity to get the desired lower 
bound part of Theorem 1. It turns out that if x points in the direction of one of 
the two vectors (ei, e-i) of the canonical orthonormal basis of Z 2 , then a judicious 
symmetry argument of Chayes, Chayes and Campanino [10] permits to obtain the 
desired lower bound. 

Proposition 6. There exists (3$ > (3 C and a positive constant c, such that for all 
(3 C < (3 < Pa, for all n > 1 and for all x = aei,i = 1,2 where a G N \ {0} the 
following holds 

,'l rrr -(8-8,)n/2 2- e -(g-feW 2 ^ ;; 

H p [a(0)a(nx)] > \ ——A— — | exp(-nr /3 (a;)), 

so that for every double sequence n f oo and (3 J. (3 C satisfying (3), we have 

l™ l -o\ iog^[a(0)a(nx)] = -4\x\. 

n,0 (p - (3 c )n 

Proof. Let n > and (3 > (3 C . By Lemma 5, we have that 

J2 A<r(0Mv)] > exp(-nT P {x)). 

yGw 13 (nx) 



On the other side, 

2 < 



3 (. V )-nTg(x)>(_0-0c)n 



< e -nr ,(x) g l{y g . fc < M^^nx) < fc + 1}| e - fe (/3-/3 c)n _ 



(/?-&)» 
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By Proposition 2, there exists /3q > (3 C such that 



(10) 



V/3 C < (3 < /3 Vz e Z 2 2\z\< 



0- & 



<Q\z\ 



Therefore for all k > 

y G w^(nx) : fc < 



< k+1 } C 



. fcn n. (fc + l)n I 
e ( — + -\x\, - — 7r -!— + 3n\x\ ] ^ . 



6 3' " 2 

From where we conclude that there exists a positive constant c = c(x) , such that 
\{y e «AnaO : fc < T -MzJI^l < fc + 1}| < cn(fc + 1). 

And hence 

^ fjP[a{0)a{y)] < cne' nT ^ ^(fc + 1) e - k(l3 '^ n 



T (3 (»)-nT (3 (x)>(/3-/3 c )ti 



fc=l 



O - p-(.P-Pc)n 
^" (W (1 _ e -( W)2 eX P(^^)- 



Thus 



£ A(0Mi,)] 



> 



T(3(v)-nT / j(x)<(/3-/3 c )n 



1 — cne 



-(/3-/3 c )n. 



2 _ e -(P-Pc)n 
(1 - e -(/3-/3c)«)2 



exp(-nr /3 (x)). 



Next, we bound the cardinality of the summation set above. Using (10), we have 

{y e w^inx) : T {y) - nTp{x) < (J3 - (3 c )n} C 
C {y e ^(ra) : |y| < (3|ar| + l/2)n} . 

Therefore |{y € w^{nx) : Tp{y) — nr^x) < ((3 — /3 c )n}\ < (6\x\ + l)n. Thus there 
exists a site y„ ji( 3 such that Tp(y n ^) — nrp(x) < {j3 — (3 c )n and 

1 - cnc- <0 -^ )n 2 ~ e _" <,3 _" ,3c '" 
(11) /W0),(«/„,)] > m + ( \~{ n W ^ )n)2 exp(-nr,(x)). 

By symmetry with respect to w^(nx), for every y„ ;/ g £ w^{nx): 

^ f3 W(y n ,/3)cr{2nx)} = ^[a(0)a(y n ^)]. 
Combining the last equality with the FKG inequality, we get that 

fjP[a(0)a(2nx)} >^[a(0)a(y n j,)}^[a(y n ^)a(2nx)} 
>^[<r(O)cr(y n ,0)} 2 . 
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Thus, we get from (11) that 



f /[a(0)a(2nx)} >^[a(O)a(y ni0 )} 2 



> 



1 C/te (l_ e -(/3-/3c)n)2 

(6\x\ + Tjn 



exp(— 2nrp(x)). 



Combining the last inequality with (6), we deduce that for any joint limit n | oo 
and (3 I @c that satisfies (3) we have that 



1 



lim , 

n,p ((3 - (3 c )n 



log[/[a{0)a{nx)] = -4\x\. 



□ 



Now we are ready to proceed to the completion of the proof of Theorem 1. 
Proof of Theorem 1 . 

Let x £ 1? be fixed. The upper bound is contained in Lemma 3. Thus, it remains 
to determine the conditions on the regime that guarantee that 



lim inf 



1 



n,p {(3 - /3 c )n 



\ogn^[a(0)a(nx)} > -4|ar|. 



From ii) of Lemma 4, There exists a positive constant K such that for all n > 

and (3 > (3 C 

A<r(0Mv)] > exp(-nr^x) - K). 

yeWP(nx) 

Fix e > 0. By Proposition 2, there exists /3(e) > (3 C such that for all (3 C < [3 < /3(e) 

7 2 M z ) 



e (4-e,4 + e). 



Therefore, for all < e < 4 and for all (3 C < f3 < (3(e): 
\y\ — \nx\ > en 



1 2e 

{i~p(y) ~ nT P {x)) > — — —nTp{x) + en((3 - (3 c )n 



4-e ' ' 16-e 

T p(y) ~ nrp(x) > e(f3 - (3 c )n. 



Using (6), we get 



A^)<v)\< E A^)<v)\ 



\y\ — \nx\ >en 



T (y)-T (nx)>e(f3-0 c )n 



< e -nMx) J- \{ y g W f \nx) : k < 



fe=i 



e(f3 - I3 c )n 



<k+l}\e- E( ^^ nk . 



By the same arguments than those used in the proof of Proposition 2, we have that 



\ y eW?(nx):k< T ^~ T ^ x) <k + l\ 
{ e((3-f3 c )n J 



< c\n(k + 1), 
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where c\ — c\{x) is a positive constant. Hence 



2 - e -e(P-p c )n 



yCWPinx) k=l 
Iv|-|nx|>en 



< c i™ e £(/3 e - e( ff- j?e) n ) 2 exp(-nrff(g)). 

Therefore 

/ 9 — p- £ (/ 3 -/ a c)n \ 

E /[^(OMJ/)] > e~* - Cl ne- £ ^fe)" _ — exp(-n^(a;)). 

|y|-|nx|<en 

And for any joint limit satisfying (3), we have 

hminf^-^-log £ A-(0)a(y)]>-A\x\. 

| w |-!nx|<en 

Since ^(m) is a cylinder of finite basis, there exists a positive constant c 2 = 02(21) 
such that 

|{y € W^(nx) : |y| — |nx| < en}| < c 2 n. 
Therefore there exists a site 2/ £i n )/3 in W {nx) such that |y| - \nx\ < en and 

So that for any joint limit satisfying (3) 

(12) liminf 1 log//[a(0)a(t/ £ ,„ ;/3 )] > — 4|ar|. 

n,/3 [p — p c )n 

On the other hand, by (6) and by the definition of W 13 (nx) we have that 

/H^tfeu)] ^ exp(-r /3 (y £i „ )/ 3)) < exp(-nr |3 (a;)). 

Combining the last inequality with (12) and using Proposition 2, we get that for 
any joint limit satisfying (3), the following holds 

(13) - Hm l0 8^(°y w)] = lim = l im -^L = 4 N . 

Next, by the FKG inequality and by translation invariance, we have that 
(14) 

^[a(0)a(nx)] > ^[a(0)a(y £tn ^)}^[a(0)a(nx - y s ^)] 

x A l/3 [ cr (0)cr(((«(a; - y E , n ,p/n)) ■ ei)ei)] n [a(O)a(((n(x - y e , n ,p/n)) ■ e 2 )e 2 )]. 
To finish the proof, we need the following lemma: 
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Lemma 7. Let e > 0. For all x £ 1? there exists 0(e,x) > (3 C such that for all 
(3 C < (3 < f3(e, x) and for all n > 1 the following holds 

{y £ W s \nx) : \yln\ - \x\ < e} C {y £ W f3 (nx) : \y/n - x\ < V '(4z + 2e)e + -}. 

n 

We postpone the proof of Lemma 7 to the end of this section and continue the 
proof of the theorem. By Propoposition 6 we have that 



fj, p [a(0)a(((n(x - y £ , n ,p/n)) ■ e l )e l )} > 

x cxp(-nT P (((x - y £ , n ,p/n) ■ e l )e t )). 
By Lemma 7, there exists /3(e,x) > f3 c such that: 

V/3 C < (3 < f3(e,x), n > 1 : ((or - y £ ,n,p/n) ■ < ^(4x + 2e)e + ^. 

Hence, for any joint limit satisfying (3), we obtain that 

liminf — 1 log//[cr(0)cr(((n(x - y e , n ,p/n)) ■ e,)e 4 )] > -4^/(4x + 2e)e. 
n,p (fj - [J c )n 

Combining the last result with (14) and (13), we obtain that for all e > 
liminf 1 logfjP[a(0)a(nx)] > -4\x\ - 8^(42; + 2e)e. 

n,/3 (p — p c )n 

Since the inequality is true for any £ > 0, we get the desired result. □ 

Proof of Lemma 7. The problem reduces to the following situation: given a circle 
C(0, \x\+e) centered at the origin and of radius \x\ +e, find a circle C(x, 5(e,x)) 
centered at x of radius 5(s, x), such that C(0, \x\ + e) r\w^(x) C C(x, 5(e, x)). Let 
us denote by a(f3, x) the angle formed by w@(x) and the line passing through x and 
perpendicular to (0,x), see the figure below. 



l-cne-C 3 "^)"/ 2 ^ 



(1- 



c)n/2)2 



(6|x| + l)n 



a(/3,x) 
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Elementary geometric considerations show that 5 can be chosen as follows 

5(x,e) = ^2 1 x | 2 sin 2 (a) +2\x\ (^e + \sm(a)\^sm 2 (a)\x\ 2 + 2e\x\ + e^j + e 2 . 

By Proposition 2, we have that a(0, x) — > 0, when J. /3 C . Thus, for every e > 
and for every x e Z 2 , there exists /3(e, x) > C such that 



V/3 C < < 0{e, x) : S(x, s) < V / 2(2e|a:| + e 2 ). 
Replacing w' 3 by W' 3 , induces an extra error of 1/n and we are done. □ 

3. ASYMPTOTICS OF CORRELATION FUNCTION AND RANDOM WALKS 

As promised in the introduction, we include in this section some explanations of 
the results of Mc Coy and Wu about asymptotics of correlations [11]. On page 305 
of [11], formula (4.38) and (4.39) give an asymptotic expansion of fi^[cr(0)a(x)], 
with x = (M, N), when < C is hold fixed and M 2 + iV 2 | oo. Actually, a closer 
look to the computations shows that (4.39) is just the expansion of the following 
double integral given by (4.22) and (4.23) of [11]. So that the results of [11] can be 
rewritten as 

[sinh- 4 (2/3) - 1]V4 y» cos(M6 1 +N6 2 ) 



3r , m [sinh- 4 (2/3) - l] 1 /* r* 



d6 1 d0. 



l-^(cos0i+cos0 2 )' 



when /3 < /3 C is fixed and \x\ 2 = M 2 + N 2 | oo. Where a/j = (1 + tanh 2 (/3)) 2 and 
7/3 = 2tanh(/?)(l — tanh 2 (0)). Instead of expanding this integral, we remark that 
this is just the generating function of a simple random walk. Indeed, let (S n ) n >i 
be a simple symmetric random walk on Z 2 starting at the origin and let me (0, 1), 
then it is known that for every x = (x\, x 2 ) € Z 2 : 

x k 1 r r Ja Ja cos(x 1 e 1 + x 2 2 ) 

\ P(S k = x)m k = —^\ \ dd x dd 2 - v , ^ttt. 

^ ^ 2 1- f(cos(0i) + cos(# 2 )) 

Hence, if we denote by £ m a geometric random variable of parameter m and inde- 
pendent of the walk then (15) can be rewritten as follows 



(16) K 'XOM-'-)] - [Slnh " 4 ^ ^ E(V m{0) (x)), 



here V m (x) = X)fc=o ^{s k =x} is the number of visits the killed random walk do to 
the site x. The surviving rate m(0) of the walk is given by the following formula 

(17) m(P)- 2 ^-- — - — 



a/3 sinh(2/3) + sirm (2/3) + sinh(2/3) 
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0.4- 



0.441 0.6 0.8 



1.2 1.4 1.6 



m(/3) as a function of f3 c < [3 < 2. 

In words, it is possible to calibrate the surviving rate of a killed random walk, 
such that the asymptotics of the two point function of the 2d-Ising model can be 
described in terms of the random walk. 

3.1. The surface tension at fixed temperature. In this section we show how 
to derive the exact formula for the surface tension from (16) using large deviation 
estimates for random walks. This will give a nice expression for involving the 
Cramer function of the random walk. 

Proposition 8. For all x e R 2 and for all (3 > f3 c , we have that 

Tp{x)=mU\{[3) 1 + 1 I{xh)), 

7>0 

where I(-) is the Cramer function of the simple random walk on Z 2 and 



=-log(m(/3)) =log 



sinh(2/3) + sinh(2/3) j 



Proof. It is sufficient to prove the result for x e Z 2 . The general case follows from 
the fact that 773 is a norm. From (1) and from (16), we know that 

Tp(x) = - lim - log E[V m r g) (nx)}, 

ri^oo n 

where m(/3) is given by (17). For every /3 > (3 C , we introduce the quantity A = 
X(/3) — — log(m(/3)). Also, for any n > 1 and x e Z 2 , we define the first time the 
random walk hits nx 

H(nx) = inf{fc > : Su = nx}. 
It is elementary to check that 



E[cM-^H{nx))] < E[V m{0) (nx)] < ^ _ E[cxp(-XH(nx))} 
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So that 

rg(x) = — lim — log.E[exp(— XH(nx))]. 

We start with a lower bound. Let us fix 7 > 0, then 

E[exp(-XH(nx))] > E[exp(-XH(nx)),H(nx) < jn] 

> cxp(—X-fn)P(3k < "fn : Sk — nx). 

Fix e > 0. We notice that if there exists k > such that Sk — nx, then for all 

y e B(nx, en) we have 

P(S k = nx)>(^j P(S k =y). 

Therefore, there exists a positive constant c such that 
\ogP(3k < -fn : Sk — nx) > — en log 4 + sup P(3k < -fn : Sk = ny) 

y£B(nx,sn) 

> — era log 4 — clogn + P(3k < -fn : Sk G B(nx, en)), 
where B(nc, en) = {y G Z 2 : \y — nx\ < en}. 

By further restricting fc to be equal to the largest integer [jn\ which is smaller 
than jn, we get that 

liminf — log£?[cxp(— XH (nx))] > — elog4 + limsup — \ogP(S\~ n \ G B(nx,en)). 

rwoo n ' n^ca n 

From Cramer's Theorem [12] for the simple random walk, we have 

lim sup — log P(S\ ~„ G B(nx,en)) > — inf 7/(21/7). 

n^oo n yeB(x,e) 

By the continuity of the rate function 

I(xi,X2)= sup ( X1X1 + x 2 X 2 - log ( -(cosh(Ai) + cosh(A 2 )) 

(Ai,A 2 )eK 2 V V 

we get that 

liminf — log E[exp(- XH (nx))] > -2s log4 - 7/(^/7) - 7A. 

n — >oo 72 

Since the result holds for all e > and for any 7 > 0, we get that 

liminf — log E{exp(— XH (nx))] > - inf (A7 + 7/(^/7)) , 

n— >oo 71 7>0 

It remains to establish the corresponding upper bound. For this, let us fix v > 
and consider 

E[cxp(— XH(nx)), H(nx) < vn\. 
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For any large integer M, we write 

\ogE[e- XH ^ nx \ H{nx) < vn] 

M-l , , 

<log V E\e- xmnx \^-vn < H(nx) < ^-vn) 

fe=0 

%^ k k k + 1 

<log > exp(— A — vn)P( — vn < H(nx) < — ——vn). 

k=0 

Let x = (xi,x 2 ) and observe that t < n(|xi| + |a; 2 |) implies that P(S t = nx) = 0. 
Thus, we can restrict the above sum to those k that are larger than or equal to 
fc = {M/v){\xi\ + |x 2 |) — 1. For M large enough, fc is strictly positive. Fix such 
an M and a k > k a > 0, we have 

k k + 1 k k + 1 

P( — vn < H(nx) < — vn) < P(3t € ( — vn, -——vn] : S t = nx). 

Notice that if there exists a t € (-^vn, ^jj-vn] such that S t — nx, then the random 
walk (S s ,s <G [jjvn, ^-vn\) is included in B(nx,vn/M). Therefore 

•P(t7^ < H(nx) < ^^-vn) <PNt G (-^vn, ^-^-vn] : S t € B(nx, ^-)) 

<P(S^ vn] e nB(x, £)), 

where \^vn\ denotes the smallest integer greater than or equal to vnk/M. By 
Cramer's Theorem, we get 

1 k k + 1 k I y \ 

liminf — logPf — vn < H(nx) < vn) < — inf — vl —, — . 

n^oo n y M ~ M '~ yeN(x.y/M)M 

Using the continuity of the rate function /, we further get that for any e > there 
exists Mo such that for any M > Mq 

t r 1 _ . k . ■ k ~\~ 1 . k I x 

hmml — logP — vn < Hlnx) < - - vn) < ——vl —. — 
n^oo n &v M ~ x ' ~ M ' ~ M \^v 

Therefore 

1 ( k k 

limsup-log-Efe - ^™^, H(nx) < vn] < - min — vX H vl 

n -oo n 0<fe<M^M M 

Since the last inequality is true for any e > 0, we obtain 

limsup-logE[e- XH( - nx \H(nx)<vn}<- inf ( 7 A + 7/(^/7)) . 

n->oo « " 0<7<j/ 

Next, we write 

E[e-^] = E[e-^\ H(nx) < H [l - > ~ 



- e. 

~M 1 ' ' 
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Using the established lower bound we have that for any v > 

1 E[e- XH ^ nx \ H{nx)> vn] . , . A . .. 

hmsup - log E[e - xmnx)] < + ( 7 A + <yI{xM) . 

The last infimum is finite and attained at a finite value 70 . Therefore, choosing a 
value of v which satisfies 

^>y(A + /(x/ 7o )), 

we obtain that 

E ! e -XH(nx) H ( nx ) > vn \ 

hm - — TT77 — d = 0- 

n^oo E[e~ XH ( nx >] 

Hence 

limsup-logS[e- Aff ^] =limsup-log^[e- Aff( " a:) , H{nx) < vn\ 

n — >oo Tl n — >oo ri 

<— inf (7A + 7/(^/7)) . 

0<7<jy 

This completes the proof. □ 

We end this section with a random walk description of the 2d-Ising Wulff crystal. 

Corollary 9. For all x — (x\,X2) € M 2 and for all [3 > [3 C , the surface tension is 
given by 

Tp{x) — xiarcsinh(s2!i) + X2arcsinh(sa;2), 
where s solves the following equation 

2^2 , . /-, I D 2^2 ■ • • 1 



1 + s 2 x{ + \ 1 + s 2 x 2 2 = sinh(2/3) 



sinh(2/3) ' 

The boundary of the Wulff crystal is given by the following set 



dWfi = ier: C(x) = - sinh(2/3) + 



sinh(2/3) 

where L{x) = E[cxp(x • Si)] = 5(cosh(:n) + cosh(x 2 )) is the Laplace transform of 
the random walk evaluated at x. 

Proof. From Proposition 8, we know that 

t p {x) = inf (A(/3)7 + 7/(x/ 7 )). 

Since the function I is even, we have 

Tp(x) = - sup( 7 A(/3) + 7/(x/7)). 

7<0 

The function / is the Legendre transform of A(x) — log(£(x)), so that 
Tp{x) = - sup sup ((x ■ y) + 7(A(/3) - A(y))) 

7<0yeR 2 

= sup $(j/,7), 
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where 

$ : fl = (-oo,0) x K 2 — ► R 

(7,y)— +7(A(/9)-A(y)). 

Hence, there exists a point (70,2/0) € such that Tp{x) — $(70,2/0) an d (70,2/0) is 
the unique solution in of the equation V<&(7, y) = 0. Note that the last equation 
is the Lagrange equation associated to the maximization of y 1— > (x ■ y) on the set 
of y e K 2 satisfying A(y) = A. Therefore 

(18) t p {x)= max (a; • y). 

j/:A(s/)=A(/3) 

Since the Wulff shape is the convex polar of the ball corresponding to the norm 
T/s(x), we get that 

Wp = {x e M 2 : A(x) < A(/3)}, 

from which we conclude the statement about the Wulff shape. The formula for the 
surface tension is obtained by solving the constrained variational problem (18). □ 

The last corollary describes the Wulff shapes of different temperatures as the level 
sets of the Laplace transform of the increments of the simple random walk, see the 
figure below 



4.5. .. 

4., 




-5 " ''-5 



Left: The graph of A(x) Right: The level sets of A(x), i.e., the Wulff shapes. 

4. The joint limit surface tension. In this section, we derive the random walk 
analogue of the surface tension near criticality. This will give us a heuristic picture 
in terms of random walks of Theorem 1 . It is an interesting question to turn this 
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heuristic description into a rigorous construction of the surface tension in a more 
general setting than the 2d-Ising model. 

The limit [3 I [3 C corresponds, in the random walk picture, to send the surviving 
probability m to 1. Thus, in order to get the analogue of the surface tension 
near criticality, we consider the asymptotics of E(exp(— \H{nx))), in the situation 
where n | oo and A = — log(m) J. simultaneously. As one may already guess, 
these asymptotics are related to the moderate deviations of the random walk. 

Proposition 10. If n | oo and A J. in such a way that 

li f m i log i = °' 

njoo Tl A A 

then 

lim — -= logi?[exp(— XH(nx))] = —2\x\. 

nToo,A|0 ny A 

Proof. From moderate deviations results [12] for the random walk we can guess the 
following lower bound: 

E[exp(-XH(nx))] > E[exp(-XH(nx)),H(nx) < 7n/VA] 

> exp(-nVA7)P(S' L7 „ /v ^ J = nx) 

where 7 is an arbitrary positive constant. We define 

— ^ n 1 



VA jny/X 
Then 

I AT r 

E[exp(-XH(nx))} > exp(-a^ 1 )P(5 LArn j = J ^ -), 

so that 

1 I N n x 

liminf — —\ogE(exp(—\H(nx))) > — 7 + 7liminf a n logP(Si N 1 = \ — - — ), 
«T°o ny A "Too n 1 a„ 7 

Then, from moderate deviations for the simple random walk, as soon as a n J. 
which is equivalent to n 2 A | 00, we obtain 

liminf — log_E(cxp(— XH {nx))) > - (7 + 7/(^/7)) , 
"Too nvA ^ ' 

where I(x/j) = \x\ 2 /-f 2 is the Cramer function of the Gaussian approximation of 
our random walk. The upper bound is proved as in Proposition 8. We get 

liminf — \ogE(exp(—XH(nx))) < — inf [ v + vl{x/v) — lim sup — i^JL^ ] 
"Too n\JX ">° y n |oo nVX J 

By imposing lim„j 00 ,\j.o log X/nVX = 0, the lower bound matches the upper bound. 
The result is obtained by computing the infimum: inf„>o I v + vl(xjv) ) = 2|x|. □ 



Let us stress out that Proposition 10 does not represent an alternative derivation 
of Theorem 1. It would have been so if we proved that in the joint limit (3) the 
asymptotic relation (15) is still valid. Actually, such a proof is possible by adding 
further restrictions to the regime (3) and it would require a non-trivial modification 
of the proof of Mc Coy and Wu [11]. This would add another layer of explicit 
computations. 
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4.1. The explicit computations formula and the method of Laplace. Until 
now we considered only rough asymptotics of V m (xi,X2) that were good enough 
to describe the surface tension. In this section, we use the method of Laplace on 
the double integral to reproduce the prefactor in the results of Mc Coy and Wu. 
This will allow us to view the Ornstcin-Zernikc behavior as a result of a Laplace 
method. First, we express V m (xi, X2) in terms of Bessel functions. For an extensive 
treatment of these special functions, we refer the reader to [27]. 

Lemma 11. For every (xi,x 2 ) e I? we have that 

f 1 T J 2 

V m (xi,x 2 ) =—. 



7rT(x 1 + l/2)T(x 2 + 1/2) 



/ du d6i / d9 2 exp(rf m (u, B 1 ,6 2 )), 
J[0,+oo) J[oM J[o,w] 



where 



TTLU 

f m (u,0- L ,6 2 ) =— (cos^i + cos 6 2 ) -u 

. f mu sin 2 Q\\ . I mu sin 2 9 2 

+ cos ip log — + sin ip log — — 

\ Acomp J \ 4sin<£ 

and r = \J x\ 2 + x 2 2 ■ 

Proof. Noticing that each walk that starts at the origin and ends at (x\, X2) has to 
contain x\ + 2k\ horizontal displacements and X2 + 2k 2 vertical displacements with 
k\ and k 2 are two positive integers, we get that 

00 

V m (x u x 2 )= ]T mXl+X2+2{kl+k2)p (X Xl+X2+2 (k 1+ k 2 )^(xuX2)) 

fci,fe 2 =0 

^ Oi +x 2 +2(fei +k 2 ))\ /mY 1+X2+2{kl+k2) 
~ . \* n A:i!fc 2 !(n + fci)!(x 2 + k 2 )\ vTJ 

Using the identity due~ ru u n = n\ r - " -1 , we obtain 

(mru/4) xi+2kl ^ (mru/4) X2+2k2 



V m {xi,x 2 ) 



Jo £± hK*i + ki)\ 



Z_, kl l {xi + kl) l fc 2 !(, 2 + fc 2 )! • 

Next, we recognize the expansion of the Bessel function of the first type [27] 

n+2k 



VneN Viet IJx) = V ,,. 1 ,., f^V 
y 1 ^ kUn + kV. \2J 

k=0 

and get 



(20) \',J-r l ...r 2 ) = r I e~™ 'I Xl (^u)I X2 (^u) du. 



Using the following representation of Bessel functions [27] 
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( /2i 

I n (an) = —= — —, / dO cxp(n(a cos6> + logsin 2 9)), 

Vtt T(i +n) Jo 



we write 



x'l 



/mru\ \Acosip) Cm r ( mu r. . 2 „ 

/ - 1 (^J = 0Fr(x 1 + i) J / o ^-pK— cos0 + cos ^ logsin 

= — ^= 7^ TT / dOcXT) 



2 

and in a similar way 



, mu uiH sin 2 I 

r — - cos + cos ip log ■ 



T fmru\ i? 

J "(—) = ^r (a2 + i/2) y " exp 



4 cos 92 



mu . musin 

r cos + sin o? log 

' 2 4 sin y 



The result follows by replacing the two last expressions for the Bessel functions into 

(20) . □ 

Lemma 11 gives a convenient way to treat the asymptotics with the method of 
Laplace. We will use the following version of this method and refer the reader to 
[14], for its derivation. 

Proposition 12. Let D C R d be an open set and consider a function f : D — > M 
that satisfies 

i) Vr > exp(r/(x)) is integrable on D. 
ii) f is twice differentiable on D. 
Hi) f reaches its global maximum on D. 

iv) The maximum of f is reached at a unique point x* € D \ 3D. 
v) Vxo e dD lim x ^ Xo f(x) = -00. 

Then one has 

(21) /cxp(r/(s))dE~f^V /2 r^oo, 
1 ' Jd 1 " \r ) y/dct(-H(x*)) 

where H(x*) is the Hessian of f at the critical point x* . 

Now we are ready to compute the asymptotics of (19), in a situation where to is 
fixed and r — > 00. 

Proposition 13. Let j3 < (3 C be fixed. For every ip e [0,7r/4], the asymptotics of 
fj,P[a(0)cr(r cos ip, r sin ip)] when r — > 00 are given by 



(sinh" 4 (2/?) - 1)V4 ^*J 1 



(1 -tanh 4 (/3)) 2 sinh(2/3) 

(22) 



4 (fiW xinhO.fo / / ; / ; \ 1 / 2 

cos 2 ^V^ + 1 + sin V^ + 1 



( ( . , 2 cos 09 . . 2sinw\\ 

exp — r cos 95 arcsinh h sines arcsinh , 

\ \ mu* mu* J J 
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where 

'((1-m 2 ) sin 2 2ip + m 2 ) L/Z + 1 



i /o x 1/2 



U = 



1 — m 



2 



and 

2 

m 



sin M2/?) + 5 i4^' 

Proof. We apply the method of Laplace given by Proposition 12 on the function 

TflU 

f m (u,6 1 ,9 2 ) =— (cos 0i + cos0 2 ) -u 

f mu sin 2 6>i \ . , / m u sin 2 2 
+ cos (p log — I + sm ip log 



4cosiyS / \ 4sin< ) 5 

defined on D = (0, +oo) x (0,tt) 2 . 

To show that f m reaches its global maximum at a unique point, we note that for 
every (61,62) € [0, n] 2 one has 

vn 

-(1 - y(cos6»i + cos0 2 )) < -(1-m) < 0. 
Thus for a; G <9D we certainly have that lim x ^ xo f m (x) = —00 because 

lim f m (u, 9i,6 2 ) = -00, 

(0i,0 2 )-(O,O) 



uniformly in u and 



lim f m (u, 61,62) = -00, 

u— »+oo 



uniformly in 0i et 2 • Hence f m reaches its maximum in D \ 3D at a critical point 
that satisfies 



-1 + f ( C OS0 1 +COS0 2 ) + ^p + ^ = 

V/ m (u,0i,0 2 ) = { -f wsin0! + 2cos(^cot0i = 

- 1« sin 2 + 2 sin <p cot 02 = 

It is easy to solve these equations and to get that it admits a unique solution 

(«*, 6l,6* 2 ) given by 



u = u 



1 + a/1- (1-m 2 ) cos 2 2ip 



+ V ~ 1-m 2 



cos 61 = ex P f -arcsinh (^) ) = J (^) + 1 - ^ 
\ \ mu* J J ]j \ mu* J mu* 

/,* / • , /2sin<p\\ /2s'm(p\ 2 2 simp 

cos0 2 = exp — arcsmh = 4/ +1 

\ \ mu* J J V \ mu* J mu* 
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The Hessian matrix H(u, #1, #2) of f m at a point (u, 9i,6 2 ) is given by 
/- cosy +/ iny -fsin^ -f sin 2 

-fsin^ 
V -f sintf 2 -m ucosd2 + 2 ^_J 

and at the critical point, the determinant of H* = H(u* , 6*, 9%) can be computed, 
we get 



2 / /n \ 2 



2 m J 2 / / 2 sin <p \ . 2 / / 2 cos (/? 



(23) det W * = --^cos^— +l + sin>y^) f 1 | < 0. 

Thus, the method of Laplace can be applied. For this, we compute the maximum 
of f m and get 

,/*/»* . , 2 cost/? . 2 sin 93 

(24) / m (u ,0 1 ,0 2 ) = — cos arcsmh sin </? arcsmh cos y> — sin ip 

rau* rau* 

The asymptotics of the integral in (19) are then obtained by replacing (23) and 
(24) into (21). We get that for each < m < 1 



(f , » , ^Try/'exp^/™^*)) 
exp(r/ m (a;))cfa; ~ — 



(25) 



[0,+oo]x[0,tt] 2 V r J ^dct(-H*) 

cxp (-r ( cos tp arcsinh 2 -^ + s i n v ^cs\nh?^-\) e r{ - cos * +sin ^) 



(^) 3/2 m 1/2 (cos 2 v vlS+i + sin2 ^V^ +1 



1/2 



when r — > 00. Finally, we treat the prefactor x^x 2 2 /nT(xi + l/2)T(x 2 + 1/2). 
From Stirling's formula, we get 

X? 1 X% 2 r re Kcos V+S in V ) 

(26) r(, 1 + i/2)r(,2 + i/2) , 2^ ' ' = 10^)1-00. 

Combining (26), (25) and (16) , we get the desired result. □ 
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